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We study the entanglement Renyi a-entropy (ERaE) as the measure of entanglement. Instead 
of a single quantity in standard entanglement quantification for a quantum state by using the 
von Neumann entropy for the well-accepted entanglement of formation (EoF), the ERaE gives a 
continuous spectrum parametrized by variable a as the entanglement measure, and it reduces to 
the standard EoF in the special case a —^ 1. The ERaE provides more information in entanglement 
quantification, and can be used such as in determining the convertibility of entangled states by local 
operations and classical communication. A series of new results are obtained: (i) we can show that 
ERaE of two states, which can be mixed or pure, may be incomparable, in contrast to the fact that 
there always exists an order for EoF of two states; (ii) similar as the case of EoF, we study in a fully 
analytical way the ERaE for arbitrary two-qubit states, the Werner states and isotropic states in 
general d-dimension; (iii) we provide a proof of the previous conjecture for the analytical functional 
form of EoF of isotropic states in arbitrary d-dimension. 

PACS numbers: 03.65.Yz, 03.67.Lx, 03.67.-a, 42.50.-p 


I. INTRODUCTION 

Entanglement is a valuable resource for quantum infor¬ 
mation processing [ij. Quantification of entanglement is 
a fundamental problem in quantum information science 
and quantum physics. Various measures of entanglement 
have been proposed such as, entanglement of formation 
ffioF), distillable entanglement and entanglement cost 
d, 3 , relative entropy of entanglement 0 , see reviews 
for more results These measures have different, yet 
closely related, physical interpretations. In general, they 
can be associated respectively with different protocols for 
quantum information processing. Several well-accepted 
measures of entanglement such as EoF and the relative 
entropy of entanglement converge to the same quantity 
for pure bipartite state, which is the von Neumann en¬ 
tropy of the reduced density operator of this bipartite 
state. For a given state, entanglement measures are not 
the same in general, nor a unique quantity even within 
one kind of measure. A class of measures may constitute 
the entanglement monotones with physical significance 
in the framework of local operations and classical com¬ 
munication (LOCC) iOl- Renyi a-entropy is a natural 
generalization of von Neumann entropy and it reduces to 
the latter when a is approaching 1. 

In this paper, we shall consider the quantity entan¬ 
glement Renyi a-entropy (ERaE) as the entanglement 
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measure, which is a generalizing of the well-known EoE. 
These entropies parameterized by a continuous variable 
a can be spectrum of entanglement monotone. Impor¬ 
tant applications are found by using Renyi a-entropy in 
describing entanglement of ground states of many-body 
systems p^ - [l9l |. which are pure states. 


In contrast to the relatively simple case of pure en¬ 
tangled states, the quantification of mixed states entan¬ 
glement is still challen^g due to the need for hard op¬ 
timization procedures [7[. However for EoE, analytical 
results are well known for some special cases, including 
arbitrary two-qubit states based on concurrence in the 
seminal work of Wootters [1^, the isotropic states [1^ 
and the Werner states (2^ in arbitrary d-dimension. In 
parallel with those three analytical results, in this paper, 
we obtain similarly analytical results of ERaE for those 
cases. A series of results are obtained, which are differ¬ 
ent from and complementary to the EoE in quantifying 
entanglement. We show that two-qubit states may be in¬ 
comparable, which is in contrast with the expectation by 
using concurrence that entanglement of two-qubit states 
can be fully quantified, which is also different from the 
case of two-qubit pure states. We can show that two 
mixed states may be incomparable. We find that the 
optimal pure states decompositions should not be the 
same, for example, in ERaE. Besides, we provide a proof 
to the conjecture about the analytical functional form of 
EoF for isotropic states. 
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II. DEFINITION 

Suppose we have a composite system with subsystems 
A and B in a pure state \il}) whose Schmidt decomposi¬ 
tion is \ij}) = bi)AB where jl is the Schmidt 

vector. For simplicity, we denote the density matrix as 
Ip = \tp){'tp\, and let Pb{A) = be the reduced 

density matrix of subsystem B (or A). The entanglement 
of pure state |^) can be quantified by the Renyi a-entropy 
of one of the reduced density operators, for example ps, 
which is defined as, 

= (l-a)”^log(trp^), (1) 

As pointed out in the introduction, TZa{ip) is reduced to 
the well known entanglement measure of von Neumann 
entropy in the a —>■ 1 limit: TZa^i{ip) = —tr(/9B logps). 
This measure of entanglement can be easily generalized 
to mixed states using the so-called convex roof construc¬ 
tion 0, [HI- For a mixed state with density matrix 
p,ERaE is defined as, 

TZa{p)= min 'S^Pk'R-aii’k)- (2) 

{pk,Pk} , 

k 

where the minimization is over all possible pure state 
ensembles {pk,'4’k} satisfying p = '^f.Pk'fpk- As in most 
cases of mixed state entanglement measure, the evalua¬ 
tion of ERaE of mixed states is much more difficult to 
carry out due to the complexity involving in the opti¬ 
mization. 

Our calculation of REaE of symmetric mixed states is 
closely based on that of two-qubit states, hence at the 
outset we introduce useful concepts and extend existing 
results concerning REaE of mixed two-qubit states to a 
broader range of a. 

It is well known that EoF corresponding to a —>■ 1 for 
ERaE depends only on concurrence which has an ana¬ 
lytic form for arbitrary two-qubit state and EoF itself can 
act as a measure of entanglement [2^. For a € (1, -l-oo), 
ERaE depends similarly only on concurrence |2l| . It is 
then possible that concurrence, in principle, might be 
the only essential measure of entanglement for two-qubit 
state even ERaE for a G (0,1) can act as entanglement 
monotones 0. We will show in this paper that this is 
not the case. We remark that ERaE satisfies monog amy 
inequality [Hill for multiqubit states when a = 2 [2l| . 
EoF or concurrence of mixed states in higher-dimensional 
system are known for classes of states with special sym- 
metry such as the Werner states [H and isotropic states 

MM- 


III. THE CRITICAL VALUE OF a FOR 
TWO-QUBIT STATE 

Suppose we have a composite system of two-qubit in a 
pure state |'0)- Given the spin flip operation on this state, 
\ip) = ay\tj;*), its concurrence can be defined as, C{tp) = 


|('0|'0)|. For notation simplicity we use C in the following 
to denote C{ijj) or C{ip) when no confusion arises. It is 
easy to see that the Schmidt coefficients X± of |^/>), i.e. 
the eigenvalues of the reduced density matrix ps, are in 
one-to -one co rrespondence with C via the relation A± = 
(1 ± Vl — C2)/2. By direct substitution of A± into the 
definition Eq. (HD of ERaE, we obtain 

7^a(V') = (l-a)-Mog(A“+A“) 

= Q(C,a), (3) 

where we introduce the function Q(C,a) for later conve¬ 
nience. 

The concurrence of arbitrary mixed state p can 
be similarly defined through the convex roof formula, 
C{p) = niin{p^^_^j^} (V'fe): where the minimiza¬ 

tion is again over all possible pure state decomposi¬ 
tion of p. An important observation made in [lol . 
Hsjl is that this measure is computable for two-qubit 
state. If we generalize the spin flip operation to any 
mixed state p hy p = {<Jy ® Gy) p* {cfy ® Gy), and let 
(Ai > A 2 > A 3 > A 4 > 0 ) denote the eigenvalues of 
pp, then the concurrence can be calculated explicitly by 
C (p) = max {Ai - A 2 - A 3 - A 4 , 0}. 

The ERaE of an arbitrary state p now becomes 

TZa{p)= min ^pkXl {C (-ipk) ,a). (4) 

iPk.Pk} , 
k 

Before we proceed, we first note some existing work on 
ERaE. It is proved in [H that ERaE satisfies monogamy 
inequality [Hj IH for multiqubit states when a = 2 . 
Also, EoF and concurrence of mixed states in higher¬ 
dimensional system are known for classes of states with 
special symmetry such as the Werner states (H and 
isotropic states [HtH]' 

Going back to two-qubit system, it is shown in (H . 
[H| that for a £ [ 1 , + 00 ) this quantity depends only on 
C(p), and in the rest of this section we investigate its 
behavior in the range of a £ (0,1). The most of the 
calculations, albeit complicating, are non-essential to out 
understanding of the results. Thus we refer the interested 
readers to appendix El and m for additional details of 
the calculation, and only discuss the indications of the 
results we obtain, which form the basis of our calculation 
of Werner states. We make use of the convexity of the 
function (C, a) which is determined from the inequality 

d^n{C,a) [<0, a£ 
dC^ |> 0 , a£ 

as well as its monotonicity with respect to C. Thus when 
a > ac = '^ 2 ^ ~ 0.82, the ERaE can be calculated 
analytically based on concurrence by 

7^Q,(p) = Q(C(p),a). ( 6 ) 

The ERaE for two-qubit states when a < ac in general 
is still a challenging problem. However, we may instead 
consider the Werner state which possesses special sym¬ 
metry [Hi- 


. 0 , 

vT-l 
2 ’ 


( 5 ) 
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IV. ERqE of WERNER STATE 

By the use of the permutation operator F = 
l*j) the Werner state pj of a bipartite sys¬ 
tem, consisting of two d-dimensional subsystems, can be 
defined as. 


^ l-i^I + F l + F’I-F 
2 d^+d 2 d^-d’ 

where the parameter F G [—1,1] specifying the state can 
be given from the relation 


tr(Fp-) = -F. (8) 

Our choice of the parameter F is different from by 
a minus sign, so that Werner states are separable for 
< 0. As we shall show later, ERaE with a = 0 of the 
qubit-qubit state is, interestingly, equal to F. 

We will make use of the result first obtained in [2^ , but 
we first introduce some notations to simplify the equa¬ 
tions. The convex hull of a function / (x) with 1) as its 
domain is defined as 


'^PkXk = x,Xk eV \ , 

k J 

(9) 

where the coefficients pk of the convex combinations sat¬ 
isfy = 1; we also need the function fw (p) = 

—tr(Fp). Thus the REoE of a Werner state is equal 
to 


CO if (x)) = inf i '^Pkfixk) 


'R-a (pD = CO (w (F, a, d )), (10) 

where the function w (F, a, d) is defined by w (F, a, d) = 
inf {Ha (tp) \fw (ip) = F, rk(V') < d}. 

To express the value of fw {"P), we write the Schmidt 
decomposition of p) as 

d d 

\P) = X! \o^, bi) = {Ua 0 Ub) ^ ^/JM K*). (11) 

i—1 i—1 

and let V = U\UB,Vij = {i\V\j), then we arrive at 

d 

fw {P) = -Y1 (12) 

*,j=i 

We first note that, since the class of Werner state is 
equivalent to isotropic states for d = 2, the calculation for 
that case can alternatively be done for the set of isotropic 
states, and the results will be identical. The function 
w (F, a) can be shown to have the explicit formula 


and quite tedious. Therefore here we only present the 
result, and the reader can consult appendix IQ for a de¬ 
tailed derivation. Substituting Ea. lfT^ into Ea. dTO . one 
concludes that 


7^a(pJ) = (14) 

[0, Fg[-1,0]. ^ ^ 

Thus as far as REaE are concerned, our result implies 
that the parameter F of Werner states is dimensionless, 
since Werner states with the same F for given a all have 
the same value of REaE regardless of the dimension of 
the Hilbert space. Also, as in the case of isotropic states, 
we obtain for the class of Werner states a relation between 
the parameter F and the REaE corresponding to a = 0 

F„=o {p7) = F. (15) 

As an example, we can compare the ERaE of a Werner 
state with F = 0.8 and a pure state with F = 0.5, as 
shown in FIG.l. We note that the pure state here is cho¬ 
sen to be a two-qubit state, in which case if we define F 
of the said state to be its concurrence, then by Eq. (|31) its 
REaE coincides in function form with w (F, a), which is 
shown explicitly in Eq. m- This will simplify the cal¬ 
culation when we perform the compaison, yet the more 
important reason for using choosing such state is as fol¬ 
lowing: as d increases the value of REaE of a pure state 
in the a —>■ 0 limit also increases, so the same behavior as 
in EIG.l will always be present as long as we set the EoF 
of that pure state to be low enough without reducing its 
d. Thus we essentially only need to verify the existence 
of a crossing for the case of Werner states and two-qubit 
pure states, which is indeed true as demonstrated in FIG. 
1: apparently, by the entanglement measure EoF corre¬ 
sponding to a = 1, the entanglement of the Werner state 
is larger than this pure state. This result seems natural 
and well-accepted. Surprisingly, when a is small and is 
approaching 0, the order of entanglement for those two 
states is reversed. One can find that the entanglement of 
the pure state is larger than that of the Werner state. 

Actually, from Eq. m we have that the entanglement 
of Werner states will always be less than that of an ar¬ 
bitrary pure entangled state in the limit of a —>■ 0. Fol¬ 
lowing the discussions above, we only need to show that 
for the case of d = 2. The optimal pure states decompo¬ 
sition for a Werner state will include a maximally entan¬ 
gled state with probability F and the identity operator, 
resulting in that ERaE equals to F, while ERaE of a 
generic entangled pure states is 1. 

In this sense, we find that ERaE of Werner states and 
pure states may be incomparable. 


H(F,a), Fg(0,1], 
0, Fg[-1,0]. 


(13) 


V. ERaE OF ISOTROPIC STATE 


Although this is crucial to our computation of REaE The class of isotropic states, specified by a parameter 
of Werner states, the derivation is essentially algebraic F G [0,1], consists of convex mixtures of a maximally 
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- Werner state 
with F=0.8 

pure state 
with F=0.5 

0.0 0.2 0.4 0.6 0.8 1.0 

a 

FIG. 1: ERqE of a Werner state and a pure state. ERqE 
of the Werner state with F = 0.8 and ERaE of a pure state 
with F = 0.5 are presented depending on parameter a as the 
x-axis. When q = 1 which is the case of the standard EoF by 
von Neumann entropy, the entanglement of Werner state is 
larger than that of the pure state. However, the order of the 
entanglement is reversed when a is close to zero. The specific 
choice of and definition of F for the pure state are discussed 
in the main text. 



entangled state and a maximally mixed state 

= FP+ + ^^^{l-P+), Fg[ 0 , 1 ], (16) 

where P+ is the projector onto the subspace spanned by 
the maximally entangled state |'I'+) = X]i=i I**)- We 
can define an isotropic state analog of fw {p) as f^i+ (p) = 
(dt+l p |4'+) = tr(P_|_p), which is the fidelity between p 
and 4'+. 

Again let us start with the formula 

TZ^{p’^°)=co{rj{F,a,d)), (17) 

where the function ij (F, a, d) is defined as ij (F, a, d) = 
inf {TZa (V’) |/>if+ (V') = iklip) < d}. Making use of the 
Schmidt decomposition in Eq. m, and define W = 
Uj^UBjWij = {i\W\j), then straightforward calculation 
yields 




1 

d 


d 

y/lpWii 

2=1 


2 


(18) 


The value of p (F, a, d) for F G [O, can then be 
easily deduced by setting = 1 , run = -^/F, which yields 
p{F,a,d) = 0. For F G ( 3 , 1 ], using the method of 
Lagrange multipliers, we derive a closed expression of 
the function p (F, a, d) as 


p{F,a,d) = - -log 7“ + (d-l)'-“(l-7)“ , (19) 

1 — a L J 


where 7 , here standing for the function 7 (F, d), is defined 
as ^{F,d) = 2 i'/F + {d — 1) {1 — F)j . In the limit 
of a —>■ 1, 77 (F, a, d) reduces to the function 


£(F,d) = F 2 ( 7 ) + (1 - 7 )log(d-1), (20) 


where H 2 (■) denotes the binary entropy function. This 
result, which we have derived here as a special case, is 
first obtained in [ 2 ^; since the calculation of p (F, a, d) is 
just a straightforward generalization, we do not go into 
the calculation details of Eq. (HU. As an application of 
the results we obtain here, we provide in the appendix [D] 
an analytical proof of the conjecture of EoF of isotropic 
states in [ 2 ^, namely 


£ (pT) 


0, 

e (F, d), 

ffl2I^(i._l) + logd, 


F G 
F G 
F G 



1 4(d-l) 
d’ IP 


4(rf-l) 
£P ■ 


( 21 ) 

Concluding this section we present an example shown 
in FIGH] obtained by numerical evaluation, of a compar- 


1.6 
1.5 
1.4 

a 

1.3 
1.2 
1.1 

0.0 0.2 0.4 0.6 0.8 1.0 

a 

FIG. 2: (color online) Gomparing the REqE with a G [0,1] 
of an isotropic state and a pure state. 


Isotropic state 
with F=0.9 d=3 

- pure state 
F=0.85 d=3 


ison between the Renyi entropy of a certain entangled 
isotropic state of F = 0.85 and a pure state chosen such 
that its EoF is smaller than that of the former. We fur¬ 
ther note that the crossing behavior is also possible be¬ 
tween two mixed states, although it does not necessarily 
have to be always present, as indicated in FIG.3: there’s 
evidently no crossing between the two isotropic states 
with d = 3, yet the two with F = 0.9, d = 2 and F = 0.7 
and d = 3 respectively have crossed REaE. 


-Isotropic state 

with F=0.9d=2 

- - - Isotropic state 
with F=0.8 d=3 

— - Isotropic state 
with F=0.7 d=3 

0.0 0.2 0.4 0.6 0.8 1.0 

a 

FIG. 3: (color online) Gomparing the REaE with a G [0,1] 
of mixed states. 



We also note that our result provides a nice corollary, 
namely the relation between the parameter F = tr (P+p) 
specifying the class of isotropic states and its Renyi en- 
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tropy with a = 0 

Fd 

7^a=o(p^) = ^i:^logd. (22) 

VI. DISCUSSIONS 

ERaE quantifies entanglement. Eor two-qubit states, 
when a > Ocj ERaE can be obtained analytically based 
on the well-known concurrence. This result implies that 
the pure states in the optimal decomposition for ERaE 
when a > Oc possess the same Schmidt vector similar 
as that for concurrence. The general analytical formula 
of ERaE even for the simplest two-qubit states is still 
a challenging problem, also the pure states decomposi¬ 
tion in general will not possess the same Schmidt vec¬ 
tor. However, ERaE for Werner states can be obtained. 
Interestingly, we notice that ERaE for the simplest two- 
qubit states may be incomparable, implying that they are 
not local convertible by LOCC. This phenomenon is pre¬ 
viously only known for higher-dimensional systems. An¬ 
alytical results of ERaE for Werner states and isotropic 
states are obtained. A series of new phenomena are found 
by using ERaE, which may stimulate more interests in 
studying quantum entanglement quantification. 
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Appendix A: Proof of the general framework of 
computing ERqE of states with local symmetry 

A general approach is derived for systematically calcu¬ 
lating entanglement monotone of states with local sym¬ 
metry in [27j |. Here we present the main steps in the 
calculations, and refer the interested reader to the paper 
mentioned above for its application in calculating entan¬ 
glement measures other than ERaE, while the latter is 
the topic of the main text. Suppose we have a set of 
states p® invariant under a group, Q, of local unitaries, 
which by definition satisfies the condition 

[p^,U®V]=Q, y(U®V)&g. (Al) 

The twirling operation is a well-established tool [1^ , 
which is defined as 

V^p= J dU{U®V)p{U®V)\ (A2) 

{U0V)eg 


Noting that twirling can be thought of as a probabilis¬ 
tic superposition of unitaries, it thus belongs to the class 
of local operation and classical communication (LOCC), 
and by definition an entanglement monotone of a generic 
state decreases or remains the same under such an oper¬ 
ation. 

By change of the integration variable, it is easy to show 
that \V^p, 17 O U] = 0 for all {U ®V) € G, which means 
that p also exhibits the symmetry described by group 
Q. Therefore, the result of twirling depends solely 
on the fidelities between the twirled state and a certain 
set of projectors, which are associated with the group of 
symmetry operations Q. This set of quantities (described 
by a vector in M) can also be used as a unique spec- 
iher of the states within the class having corresponding 
symmetry. For isotropic and Werner states, the set of 
projectors can be simplified to only one operator. Since 
these are the only sets of states we consider in this paper, 
we can restrict our discussion to this special case, where 
we define a function fg (p) with the property that 

■P‘‘p = 4.„,,y (A3) 

and the function fg (p) is a one-to-one mapping from the 
class of symmetrical states to the set of real numbers (or 
in general, real vectors). 

Now we move on to the calculation of an entanglement 
monotone: suppose X (?/>) is defined on the set of pure 
states, where stand for the correspondir^ density oper¬ 
ator, and we use convex roof construction [7| to generalize 
this measure to mixed state p by 


A (p) = inf I ^PkX {fjk) 

[ k 


'^Pkf’k = , 


(A4) 


Then by defintion of the infimum, we have 


'^PkX {ifk) 

1 k 

II 

1 '^PkX (V'fe) 

1 k 



A” (p^) > inf i ^pkX {-ipk) 


'^Pkfg {'(fk) = ^ ( • 

k } 

(A5) 


Let be a pure state ensemble which is the 

optimal decomposition achieving the infimum in the right 
hand side, namely 




^p® A = inf i ^pfcA ifjk) 


'^Pkfg = F 

k 

(A7) 
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It follows that the twirling of the new state 
Y.pTi’T is just 


k 


p' = 




( 0 )„,,( 0 ) 
k 


PF=fg(p') - P 


F=Epi”Vs(v'r) 


= p%, 


(A8) 


and from the fact that twirling opration is a LOCC op¬ 
eration, we obtain the inequality 

X {p') > X {V^p') = X {p%) . (A9) 

On the other hand, by the property of the convex roof 
construction in Ea. (IA4ll . we have 


Y.X'x (d“') > ^ feprtf') = ^ (p ). (Alo) 

k \ k / 

which combined with Eqs. (IA5l) . (IAtI) and (IA9I) provides 
a tight bound for X {p%) for explicitly calculation 

x{p%) = x{p') = Y,p^^^x[^f'>) 

k 

^Pkfoi^k) = f \ . (All) 


= inf < '^PkX {ipk) 


Now we’ve derived a directly calculable formula of 
X {p%), which can be computed via a two-step proce¬ 
dure. The first step involves an optimization over pure 
states with the same value of the function fg (p). 


x{F)=M{X{i;)\fg{^)=F}. (A12) 

While this function is evaluted in for entanglement of 
formation of both Werner and isotropic states, we note 
here that it is not immediately clear how it should be 
computed for any particular class of states or type of 
entanglement monotone. The calculation still needs to 
be done explicitly for specific combinations, given the 
general framework. 

The next step is theoretically more easy and com¬ 
putable in general, but in practice does not always yield 
elegant or explicit results without numerical calculation. 
We define the convex hull of a function f (x) defined on 
V to be 


co{f{x)) = 

'^PkXk = x,'^pk = l,Xk eV L 
k k } 

(A13) 

and the entanglement monotone X (p%) of state can 
be shown to be 


inf E Pkf (Xk) 


X{p%)=co{x{F)). (A14) 


Appendix B: Proof of convexity property of the 
function (C, a) 


The function fl (C, a) is defined as 

ll(C,a) = (l-a)-Mog(A“+A“), (Bl) 


where X± = (1 ± Vl — C'^)l‘2.. Let us present the first 
derivative of 14 (C, a) with respect to C, dfl/dC, in terms 
of the Schmidt coefficients A± by using Eq. (EU, 


d^jC^ (X) _cr_ r\a— 1 dA-i- \g—i dA— , 

dC “ (l-a)(A^-f A“)^ + dC dC 

(B2) 

The derivatives of A± are. 


dA± C 

”d^ “ 2(1-2A±)’ 
d^Ai _ 2 dA± 

dC 2 “ C 2 (l- 2 A±)”d^' 


(B3a) 

(B3b) 


Also we know, dA+/dC + dA_/dC = 0, and we intro¬ 
duce the notations, Di = |dA+/dC| = Cl2\/1 — and 
X = A_/A+. Substitute Ea. (IB3l) into the derivative in 
Ea. (IB2l) . we arrive at 


dn (C,a) 


aX^-^Di 

(a - 1) (A^ -f A“) 


(l-a;“-i) > 0 . 


(B4) 


Namely 11 (C, a) is a monotonically increasing function 
with respect to C. 

The evaluation of the second derivative of 11 (C, a) is 
more complicated, 


9^11 (C,a) _ aX^^-^Dl 

(l-a)(A^-f A“)^ 


(B5) 


where 


K 

gix,a) 


1 - - (2a - 1) (1 - x) x^-\ (B7) 


and X € (0,1] and a S [0,1). Observing that 

{k > {I-x°‘-^Y >0, ae[ 0 ,i], 

[K < {1 - ae(il), 


(B8) 


we thus conclude that the function ll(C,a) is concave 
with respect to C for a G [O, 4], 

8‘n{C,a) ^ ^ 

This result is actually opposite of the convexity of the 
function H (C,a) for a G (1, -foo). We thus give a nega¬ 
tive answer for the holding of the relation 


It is easy to see that this step is computable in general. 


7eg(p) = ll(C(p),a), (BIO) 
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FIG. 4: (color online) By using condition ^ = 0, we 

can find the dependence of Co, which satisfies this equation, 
on a. 


for a G [ 0 , i]. 

Next, we consider the region a G l)- For a fixed 
a, the second derivative of il{C, a) may have a zero cor¬ 
responding to Co in the interval Cq G [0,1]- Numerical 
calculation shows that the value of Co increases monoton- 
ically with respect to a, see FIG. 1 for the dependence 
of Co on a. So there may exist a critical value of a cor¬ 
responding to C = 1 such that the second derivative of 
n (C, a) is zero. Such a critical value ac does exist, such 
that the simplification, meaning the holding of Eq. (| 6 ]) , is 
still valid for any a larger than this value. In fact, it is 
not difficult to obtain the value of ac analytically. One 
simply considers the limit C —^ 1 and the requirement 
that. 


lim 

C-J-l 


a2f](C,a) 

5C2 


> 0 , 


(Bll) 


Appendix C: Derivation of the closed form of the 
function LO {F, a, d) 


Here we give a detailed proof of the formula 


uj {F, a, d) 


n{F,a), F> 0, 

0, F < 0. 


(Cl) 


The value of uj{F,a,d) for F G [—1,0] can be easily 
obtained by setting fii = l,Uii = Vf in the explicit 
expression of the function 


d 

fw (V') = - X! (C 2 ) 

Gi=i 


giving UJ (F, a, d) = 0, which reproduces the sep¬ 

arability condition of Werner state. While the 
optimization in the definition of uj (F, a, d) = 
inf {72.Q, {ip) \ fw {ip) = F, rk (■)/>) < d} for general d 
cannot be carried out in a straightforward fashion, it 
is rather simple for two-qubit state. Setting d = 2 in 
Eg. (lea. we have that for F > 0, 

C {iP) = 2Vmh^ > F, (C3) 

where equality holds if and only if p,i _2 = A±. Since 
the function H (F, a) is monotonically increasing with re¬ 
spect to F, we thus obtain the inequality ft {C {ip) , a) = 
T^a ('*/’) ^ (F, a). It then follows that. 


UJ (F, a, d) 


n{F,a), F > 0, 
0, F < 0. 


(C4) 


When seen as a real vector, fl = [ni] belongs to the con¬ 
vex set of real, at most d-dimensional vectors such that 
J2i=i di = l(di > 0), which we denote by Kd- It follows 
that fw ('0) can also be seen as a function fw (ft, V) of 
a real vector fi and a unitary matrix V. For F G (0,1], 
we can calculate the quantity defined by 


which is equivalent to lim K < 0. Referring to the defi- 

X —>-1 

nitions of K and g{x, a) in Eas. (IB 6 l) and (IB7I) . we derive 
the following inequality, 

[3 (a - 1) + (2a - 1) a] < 0. (B12) 

o 

The value of ac can be calculated by considering the con¬ 
dition for equality in the above expression, which gives 
us, 

a, = 0.82. (B13) 

This solution is consistent with our numerical result pre¬ 
sented in FIG. I. 


T (F, a, d) = inf d ^ | , 

(C5) 

from which the function uj (F, a, d) can be obtained from 
UJ (F, a, d) = log [r (F, a, d)], by the monotonicity of the 
logarithm function. Noting that the minimization in Eq. 
(IG5I) . with d larger than 2, always covers all possible 
combinations of } with d = 2 , we have 

T (F, a, d) < r (F, a, 2). (G 6 ) 

Physically this implies that the amount of entanglement 
contained in Werner states is upper bounded by that of 
a maximally entangled pair of qubits (for example the 
Bell states). A consequence of this constraint is that 
when equality in the infimum in Eq. (IG5I) is achieved, 
the largest /i^, which we denote by /Xmax, must be bigger 









or equal to In fact, if /imax is smaller than i, we will 
derive a contradiction of Eq. (IC61) 

Mmax < i ^ T(E,a,2) < T(l,a,2) < (C7) 

i 

Thus we have a constraint on ^max that reads 

Mmax > (C8) 

Our next step is to seek an upper bound of ^maxj which 
combined with Eq. (108^ will further constrain the pos¬ 
sible values that uj {F, a, d) may take. First we derive 

I 


some inequalities that the set of variables must 

satisfy. Noting that X]f=i have 

d 

1 = 

z=l i,j 

~(^im "k ^ ^ ^ = 1, (C9) 

i/im 

where is one of the indices such that = Umax- 
With the use of Cauchy-Schwarz inequality, it follows 
that 


d 


i¥=3 


2 

( ^ 

\ 

( - 

] 



\ 

< 

Y 


E l^'C, 

2 

= 

E 




/ 

i¥^3 




/ 


2 


j¥=im 





(CIO) 


(Cll) 


Thus the following inequalities can be easily derived 




< Y 

i¥=3 


3¥=ira 


< v/st. 


where for later convenience we set 

s=Y d-im\voiS, 

3^im 

3=i^im 


(C12) 


(C13) 


(CM) 

(C15) 


Making use of the condition in the infimum Eq. (IC5I) , 
namely fw (m, V) = Yi=i = 

—F, we have for F G [0,1] that 


<2 


d 

'y \J Mim Mi j 

3¥=im 


d 


i^3 


— 1 

<1 - 


= 1 - 


- Y + 2 

3¥^im 


d 

3¥^im 


- i: 






H- 




'\ 


t ( Y 

j^'^m 


(C17) 


To further evaluate the expression in Ea. (IC17E we note 
that this is essentially a quadratic function, with the ar¬ 
gument satisfying the inequality 

* = Y ^ Y = 

3¥=im 


d 

F + Y^^^\'^Y 
2=1 


d 

i¥^3 


(C16) 



d 


d 

<2 


+ 

Y ^/^^i^^3'^3iY 


3¥^im 


i¥^3 


and we consider two possible cases of \/t to derive a uni¬ 
form bound. First, if 

^ ^ I — I — ^ J (^ 79 ) 

3¥^im 

I 1 2 

then we have rJ)„j | < 1 - < Mim “ 
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such that 

d 

F + < 2\//iniax (1 - Mmax)- (C20) 

If instead the variables {fdi^Vij} satisfy 


d 

Mim kjim 1 ^ = k*mim 1 ^ < 1 “ I (^ 21 ) 

j¥^im 

then we obtain, similarly, 

<2(1-Ai*™), (C22) 

d 

^ ^ l^^mi I — ^'^rn l^'tm I ^ f 

j^im 

(C23) 

d 

=^F + ^ ^ IJ^i\Vii\ < 1 — I < 2 (1 — , 

i=l 

(C24) 

d 

^ < 2'\/^niax (1 - Mmax). (C25) 

Combining the two cases, and noting that = Umax, 
we derive a constraint of /imax in terms of F 


d 

F < F + y^^fii\vu\'^ < 2y/imax (1 - /^max), (C26) 

2=1 


Let n (F) = (which for simplicity we will denote 

by /r when no confusion arises), and Eqs. (IC 8 I) and (IC26I) 
give ^max < ^^{F)■ Making use of the Schur concavity of 
the function r (F, a, d), we have 


r (F, a,d) = + (1 - Mmax)" > r (F,a, 2). 

(C27) 

On the other hand we have already obtained the inequal¬ 
ity (IC 6 I) . and together the two bounds yield 


T (F, a,d) = Y^ /r“ > r (F, a, 2) > r (F, a, d) , (C28) 

2 

indicatiirg that the bound given in Ea. (IC 6 l) is indeed 
tight. We have now succeeded in calculating the func¬ 
tions T (F, a, d) and lu (F, a, d) 


T{F,a,d) =T{F,a,2) (C29) 


ui (F, a, d) 


n{F,a) F > 0 
0 F < 0 


(C30) 


Appendix D: Proof of the conjecture on EoF of 
isotropic state 

We first state the rigorous result on entanglement of 
formation of isotropic state Pp ° in d dimensional Hilbert 


space, which we denote by £ (pf°) 

£{pT)=co{e{F,d)) (Dl) 

where co (•) denotes the convex hull of a function, and 
the function e (F, d) is found to be 


e (F, d) = H 2 (7 {F, dj) + (1 - 7 {F, d)) log {d - 1) (D2) 
7 (F,d) = I(v^+y(d-l)(l-F))' (D3) 


and H 2 (x) is the binary entropy function. 

In Vollbrecht and Terhal’s work [ 2 ^, a closed expres¬ 
sion of EoE of isotropic states is conjectured, whose va¬ 
lidity the authors argue can always be verified for any 
given d, by directly plotting the function e (F, d) and ex¬ 
amining its behavior. The rigorous proof of d = 3 case is 
provided. It is now generally accepted that the conjec¬ 
ture is true for arbitrary d, it seems that a proof is still 
necessary. Alternatively, here we seek to prove this con¬ 
jecture in an analytical fashion, without any presumption 
about the value of d. 

The proof consists of two steps: we first prove a gen¬ 
eral statement about £ (F, d), when treating F as the 
argument and d as a parameter, and the next step sim¬ 
ply involves the verification of a criterion of the point 
F = First we wish to show that, the concavity of 

e (F, d) with respect to F on the interval F S [1,0] is in 
general analogous to that of the special case with d = 3, 
i.e. first concave upwards and then concave downwards. 
We will prove this statement by directly evaluating the 
value of second derivative of e (F, d) with respect to F. 
To keep the equations as simple as possible, we use e and 
7 to denote respectively the functions 


e (F, d) = F 2 ( 7 ) + (1 - 7 ) log (d - 1), 
7 (F,d) = I[VF+v'(d-l)(l-F)'' 


A bit of algebra gives us 


dj ^ ^ 7(1 - 7 ) 

dF ^F(l - F)’ 

d^7 y/d — 1 ( 

dF^ ~ 2d7(l-7) v'^(l-F) \dF) 


(D4) 

(D5) 


(D 6 ) 

(D7) 


Substituting the derivatives of 7 with respect to F into 
the relation 


d'^e d'^e /' d'y\^ d ^7 de 

^ '^ydFj 


(D 8 ) 


we obtain 


d^e _ ^/d^ 7(d- 1) _ 2d^/F(l-F) 

2d[F(l-F)]5 1-7 Vd- 1 

(D9) 

Now we only need to examine the sign of the term in 
the square bracket of the above equation, which is the 
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o2 

same as that of In a change of variable, we let 

X = , so that X S [0,1] decreases monotonically 

with respect to F in the interval F G [^) l] • Let the term 
in the square bracket of Eq. (ID91) be a function / (x) of 
X, we have 




dx 


+ (d — 1) x^ ’ 


(DIO) 


whose behavior at the two ends of the interval x G [0,1] 
is easily found to be 


F 


1 


, X —>■ 1, 


F —>■ 1, X —^ 0, 


52 £ 

9F2 ^ 

(Dll) 

d^s 

5F2 ^ 

(D12) 


Differentiate / (x) with respect to x, we obtain 
{d — 1) (d — 2) dxg (x) 


/' (x) = 


(1 — x) {dx — X + 1) (1 + (d — 1) x^)" 


, (D13) 


whose denominator is non-negative, and the function 
g (x) is equal to 


9{x) 



d? 

{d-l){d-2f' 


(D14) 


where the function h (x) is defined as 

/i(x) = Inx-I-b (D18) 

For d > 2, we have h' (x) = — < 0, and it follows 

that 


f{x+)<Q. (D19) 

Combining this with the special case / (0) = 0, we have 
/(x) <0 and < 0 for x G (0,x+). On the other 
hand, the function / (x) increases monotonically for x G 
[x+, 1), and it can be easily shown that / (1) —>■ -boo as 
X —>■ 1, which ensures that the function / (x) has only one 
zero in the interval x G (0,1). Also taking into account 
Eqs. (IDTTll and dnn, which gives the limit of as 
F —>• i and 1, we conclude that the second derivative 
has only one zero on the interval F G [ 3 , 1 ]- 
Now that we have proved that the function e (F, d) on 
the interval F G [1,0] is first concave upwards and then 
concave downwards as F increases, the next step is sim¬ 
ply to find a line both tangent to e (F, d) and passes the 
point F = 1; in other words, we want to find a solution 
Fo G (^, 1 ) to the following equation 


Noting that g (x) is a quadratic function of x, it has two 
zeros which can be easily found to be 


—2\/ d — 1 i d 
(d-2) Vd^' 


(D15) 


The zero x_ < 0 is not in the interval x G [0,1] of our 
interest, while it can be shown that x+ G (0,1) if d > 
2. Furthermore, by evaluating the value of g (x) in Eq. 
(ID14I) to determine the sign of /' (x), we have 


/' (x) 



X G (0, x+) 

X G [x+,1) ■ 


(D16) 


Thus we conclude that the minimum of the function 
/ (x) on the interval x G (0,1) is / (x+) 

/ (x-^) = In Vd - 1 - = h(yd-l^ , (D17) 


logd-e=(l-F)—. (D20) 

Taking the convexity of the function e (F, d) into ac¬ 
count, one can deduce that such solution does exist and 
is unique. In addition, direct calculation shows that 
the conjecture F = is indeed the solution to 

Eq. (ID20I) . Since e{F,d) can be easily shown to be 
monotonously increasing with respect to F, we now have 
adequate information to determine explicitly the value of 

f (pr) = (pr), 



fo 

F G 


£{pT) = \ 

e{F,d) 

FG 

1 4(d-l) 
d’ £p 

1 

£ll2S^(F-l) + logd 

F G 

4(d-l) , 

^ 


(D21) 
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